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Abstract. Among those objects formed out of collisional debris during galaxy mergers, the prominent gaseous 
accumulations observed near the tip of some long tidal tails are the most likely to survive long enough to form 
genuine recycled galaxies. Using simple numerical models, Bournaud, Due & Masset (2003) claimed that tidal 
objects as massive as 10 9 Mq could only form, in these simulations, within extended dark matter (DM) haloes. 
We present here a new set of simulations of galaxy collisions to further investigate the structure of tidal tails. First 
of all, we checked that massive objects are still produced in full N-body codes that include feedback and a large 
number of particles. Using a simpler N-body code with rigid haloes, we noticed that dissipation and self-gravity 
in the tails, although important, are not the key factors. Exploiting toy models, we found that, for truncated 
DM haloes, material is stretched along the tail, while, within extended haloes, the tidal field can efficiently carry 
away from the disk a large fraction of the gas, while maintaining its surface density to a high value. This creates 
a density enhancement near the tip of the tail. Only later-on, self-gravity takes over; the gas clouds collapse and 
start forming stars. Thus, such objects were fundamentally formed following a kinematical process, according to 
a top-down scenario, contrary to the less massive Super Star Clusters that are also present around mergers. This 
conclusion leads us to introduce a restrictive definition for Tidal Dwarf Galaxies (TDGs) and their progenitors, 
considering only the most massive ones, initially mostly made of gas, that were able to pile up in the tidal tails. 
More simulations will be necessary to precisely determine the fate of these proto-TDGs and estimate their number. 

Key words. Galaxies: formation - Galaxies: interactions - Galaxies: dwarf - Galaxies: halos - (Cosmology:) dark 
matter 

1. Introduction (e.g., Ide Griis et all 120031: iTran et al.l[2003[) . Their blue 

luminosity, inferred mass and size are such that some 

Numerous papers have reported the formation of presum- of them may evolye into Globular Cluster s with typical 

ably bound stellar structures around interacting systems. masses of 10 6 M Q (e.g., ISchweizer et all fl99fih . Some 

Their blue colors and high metallicity indicate that they 1Q 7 Mq m fragments disco vered around merge rs could 

were formed rather recently out of material processed in be their gaseous progenitors fenglish et allEooi . 

the disk of their parent galaxies. The census of young _ Gmnt HU Complexes (GHCs), i.e. star-forming regions 

objects in mergers is still on going and their taxonomy with i umirLOS ities exceeding the Giant HII Regions found 

not yet completed. m i so i a t ec [ spiral disks, were discovered in long optical 

- Super Star Clusters (SSCs), with typical masses of tidal tails ( e g ^ Weilbacher, Due & Fritze-v. Alvensleben 

10 5 M Q , appear on optical images as compact stellar con- 2003; L6 p e z- Sanchez, Esteban & Rodriguez 2004). Further 

densations. They are either located in the central regions away Intergalactic Emission Line Regions, characterized 

ofmergers (e.g |Holtzman et al.| [1993 |Whitmore et al.| by a very low unc i er l y ing old stellar content, and some- 



1999t|Zepf et al.||1999|), in tidal bridges or along extended times by their compact aspect ( thc S0 -called EL-Dots) 

tidal tails (e.g., Gallagher et al. 2001; Knierman et al. wer e detected in several clusters of galaxies. Their optical 

2003; Saviane, Hibbard fc Rich 2004) and may be dust spe ctra are typical of star-forming HII regions and their 

enshrouded (|Gilbert et alj|2000|). Young Massive Clusters rather high oxygen abundances indicate that they are 

(YMCs) were identified in high-resolution HST images formed of pre -enriched gas probably stripped from parent 

colliding galaxies fe.g-. lGerhard et afll2002i [Cortese et alJ 
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2004- iRvan- Weber et alJ 12004 iMendes de Oliveira et all 



2004). 



- Finally, the formation in tidal debris of Tidal Dwarf 
Galaxies (TDGs), i.e. objects with apparent masses and 
sizes of dwarf galaxies, has been reported in several 
studies of interacting; sys tems fe.g.. lDuc fc Mirabe]| l"l999: 
iHibbard fc Barnesl I2004L and references thereein). They 
contain large quantities of gas in atomic, molecular and 
ionized form and have l uminous masses of typically 
10 9 M fe. g .. lRrainee raTl l200lh . Because they are most 
often observed at or near the tip of long optical tidal 
tails, the very existence of such massive objects has been 
challenged. Indeed an apparent accumulation of tidal 
material could be artificial. In 3-D space, tidal tails are 
curved. Seen edge-on, they appear as linear structures 
and may show at their tip fake mass concentrations 
caused by materia l proje c ted along the line of sight 
l|Hibbard fc Barnesl 120041: iMihoj l2004|) . Kinematical 
studies of tidal tails help in identifying such projection 
effects l)Bournaud et alJ 12004) . The spatial and velocity 
coincidence of the different phases of the tidal gas towards 
a TDG candidate adds ci rcumstantial evidence as to its 
reality l|Braine et al J I2001 1 . 

These objects around mergers all contain stars with an 
age, derived from their colors or spectroscopic properties, 
which is comparable to or smaller than the dynamical age 
of the collision between the parent galaxies. They were 
hence most probably formed during the merging process. 
This does not mean however that a single mechanism is 
at their origin. Indeed, the heterogeneity of their prop- 
erties (gaseous content, stellar population, location, and 
above all their luminosity/mass) rather suggests that sev- 
eral physical processes take place during the collision. The 
initial conditions of the encounter may also matter to ac- 
count, for instance, for the formation of the most massive 
tidal objects, exceeding 10 9 M©, i.e. those often referred 
as Tidal Dwarf Galaxies. 

Numerical simulations of collisions are particularly 
helpful to identify the mechanism responsible for the 
formation of young objects in tidal debris. Soon af- 
ter the study by Mirabel, Dottori & Lutz (1992) of a 
Tidal Dwarf Galaxy candidate in the Antennae system, 
iBarnes fc Hernauistl l|l992|) and Elmegreen, Kaufman & 
Thomasson (1993) published numerical models exhibit- 
ing bound condensa ti ons a long tidal tails. While in the 
IBarnes fc Hernauistl l|l992l ) simulations, the condensa- 
tions formed from gra vitational instabil i ties in the stel- 
lar component, those of lElmegreen et al 1 lll993l) were pro- 
duced directly from gas clouds, the velocity dispersion of 
which had increased due to the collision. Further simula- 
tions, made with a higher resolution and/or more realistic 
initial conditions, especially regarding the respective dis- 
tribution of the gaseous and stellar components, could pro- 
duce along the tails objects lik ely to be the progenitors of 
SSCs or globular clusters (e.g.. IBarnes fc Hernauisti ri996: 
lKroupall99llBekki et all2002h . Their masses of typically 
10 s M are however one order of magnitude lower than 



that measured in the Tidal Dwarf Galaxies located at the 
tip of tidal tails. 

In a few older simulations, however, the formation of 
massive a ccumula ti ons o f tidal material were noted (See 
Fig. 3 in IBarnes! l|l992|) or Fig. 3b in lElmegreen et alJ 
l|l993|) 1 but were not much commented as they were 
quickly disrupted or went out of the volume studied in the 
model. lElmegreen et al.1 l|l993h actually concluded that 
massive galaxy interactions could lead to the ejection of 
the most outer HI disk, forming an "extended 10 9 M© pool 
at the end of the tidal tail" . Recently, Bournaud, Due fc 
Masset (2003) (hereafter BDM03) using a new set of sim- 
ulations involving simple co-planar collisions showed that 
such large end-of-tail blobs could condense and become 
gravitationally bound. Varying the initial size of the dark 
matter (DM) haloes in the parent galaxies, they found 
that objects as massive as 10 9 M©, i.e. similar to TDGs, 
could only be produced when the haloes were extended 
enough - at least ten times the optical radius. Most often, 
dark matter haloes are truncated in numerical simulations 
although cosmological models predict their large size. This 
trick saves a great number of DM particles and is justified 
as long as the central regions of colliding galaxies are stud- 
ied. BDM03 actually found that the extent of the DM halo 
influences the internal structure of the tidal tails, i. e. the 
distribution of matter along them, while previous works 
had more focussed on the length and curvature of these 
tails (e.g., Mihos, Dubinski & Hernquist 1998; Springcl & 
White 1999; Dubinski, Mihos & Hernquist 1999). 

The simulations presented by BDM03 had a number 
of limits. The DM and stellar components were computed 
with a 3-D FFT code at a resolution of only 5 kpc while 
the gas was simulated in 2-D (i.e. in a co-planar en- 
counter) with a resolution of 150 pc high enough to re- 
solve the gaseous tidal tails. Thus, the possible role of 
stellar sub-structures (which could not form in the low- 
resolution code) in building TDGs was ignored. Beside, 
the gas/DM halo interaction was not taken into account. 
In this paper, we check the results of BDM03 in full 3- 
D N-body simulations that have, beside, a much higher 
resolution in both the DM, stellar and gas components 
(Sect. I2J . Having confirmed with this new code the early 
results of BDM03 and assessed the fundamental role of 
kinematical processes in structuring tidal tails, we use sim- 
pler simulations and qualitative arguments to clarify the 
effect of tidal forces within extended DM haloes (Sect. OH). 
Finally, in Sect. 01 we revisit the definition of Tidal Dwarf 
Galaxies, arguing that, being formed according to a top- 
down scenario, they have a different physical origin than 
the other tidal objects. 

This paper only discusses the physics of the formation 
of massive TDGs. It does not address their abundance nor 
their evolution and survival time. Answering these issues 
of cosmological importance requires to carry out a large 
number of simulations with a comprehensive set of initial 
conditions. This work is in progress and will be presented 
elsewhere. 
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Run 


N a 


M" 


m g c 


m/ 


V e 


b> 


h 9 


t 2 h 




A 


1024 a 


1.2 


8% 


10% 


150 


90 


15 


20 


27 


B 


512 3 


1.7 


5% 


12% 


150 


50 


45 


20 


40 


C 


512 3 


1.0 


8% 


8% 


170 


75 


15 


15 


30 


D 


512 3 


1.3 


9% 


9% 


120 


120 


25 


25 


35 



a Spatial resolution: N is the number of cells on the Cartesian grid. See text for the 
associated number of particles and softening length. 
Stellar mass ratio. 

c Gas mass fraction in the most massive galaxy; m g is the gas to disk+bulge mass ratio. 
d Gas mass fraction in the less massive galaxy 

e Relative velocity of the two galaxies computed at an infinite distance in km s _1 . 
' Impact parameter in kpc. 

9 Inclination of the most massive galactic disk with respect to the orbital plane. 
h Inclination of the less massive galactic disk with respect to the orbital plane. 
8 Angle between the two galactic disks. 

Table 1. Run parameters 



2. Formation of TDGs in full 3-D N-body 
simulations 

We have carried out a series of full 3-D N-body simula- 
tions of the merger of two spiral galaxies with mass ratios 
between 1:1 and 1.7:1 and different orbital parameters. 
Their resolution, which was much higher than in BDM03, 
allowed to precisely study the internal structure of the 
tidal tails. 



2.1. Code description 
2.1.1. Physical model 

We have used t he N-body FFT code of 
iBournaud fc Combesl ||2003(k In the main simulation 
presented in detail hereafter, the number of particles for 
the most massive galaxy was 2 x 10 6 for dark matter, 
10 6 for stars, and 10 6 for gas. For the other galaxy, 
the numbers were scaled to its mass. The gravitational 
potential was computed on a Cartesian grid of size 1024 3 , 
with a softening length of 390 pc. We have also run three 
other simulations with a four times smaller number of 
particles meshed on a 512 3 grid, and a softening length of 
780 pc. 

The dissipative nature of the I SM has been modeled 
throug h the sticky-particles code of IBournaud fc Combesl 
(2002J), with elasticity parameters p\ — (3 r — 0.8 for cloud- 
cloud collisions. Star formation was described by a gener- 
alized Schmidt law l|Schmidtlll959j) . and we employed a 
time-dependent stellar mass-loss model inspired of that 
of Jungwiert, Combe s fc Pa lous (2001) and described in 
IBournaud fc Combesl (2002). We then computed the star 
formation and mass-loss rates in each cell of the Cartesian 
grid, according to the values of these two rates, convert- 
ing a fraction of gaseous particles into stellar ones and vice 
versa. 



2.1.2. Mass distribution of galaxies 

Each galaxy was made-up of a spherical dark halo, a stel- 
lar and gaseous disk, and a spherical bulge. For the most 
massive galaxy, with a visible mass of 2 x 10 11 Mq, the 
stellar disk was a Toomre disk with a radius of 15 kpc, 
a radial scale-l ength of 5 kpc and the sam e vertical dis- 
tribution as in IBournaud fc Combesl l|2002j) . The gaseous 
disk is a Toomre disk of truncation radius 37.5 kpc and a 
radial scale-length of 15 kpc. Such a gaseous disk which 
is 2.5 times more extended than the s tellar one is repre- 
senta tive of most spiral galaxies (e.s.. iR.oberts fc Havnesl 
1994). The bulge was modeled by a Plummer sphere of 
scale- length 1 kpc and truncation radius 3 kpc. The dark 
halo was a softened isothermal sphere truncated at radius 
150 kpc (10 times the stellar disk radius). Its density pro- 
file is given by: 



We used a core radius r c = 4.5 kpc. For the other galaxy, 
the radial distribution of matter was scaled by the square 
root of its stellar mass. In all simulations, we chose a bulge 
to disk mass ratio of 0.3, and a dark-to- visible mass ratio 
inside the stellar radius of 0.6. The mass fraction of gas 
inside the disk has been varied (see Table 0). 

2.1.3. Galactic encounter and orbital parameters 

Tabled lists the values of the orbital parameters used in 
our simulations. The ratio M of the stellar masses between 
the two galaxies roughly corresponds to the ratio of the 
dark halo masses and hence of the total masses. We chose 
orbital parameters that favor the formation of a long tidal 
tail for at least one of the two galaxies; in particular we 
have only considered prograde orbits. 
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Fig. 1. Evolution of the gas distribution for the full N-body simulations of run A. The system is seen face-on. The 
gas is displayed with a logarithm intensity scale. 



2.2. Results 

Fig. H presents the evolution of the gas distribution for 
simulation A. The stellar, gaseous and dark matter compo- 
nents are shown at key moments of the merging sequence 
in Fig. Hand Fig. [3 

In run A, two major tidal tails of ~ 100 kpc developed 
soon after the periapse (at t =410 Myr). At t =470 Myr, 
massive gas clouds have piled up at the tip of each tidal 
tail: they are the progenitors of the TDGs. We find that 
70% of the gas along the tails is located in the 20 last 
kpc, while the tails are about 100 kpc long . The total 
masses of the two TDG progenitors at t =470 Myr, com- 
puted in spheres of diameter 20 kpc centered on the TDG 
luminosity peak, are resp. 2.9 xlO 9 M Q and 2.3 xlO 9 M Q 
(see Table 0). At this stage, most of the mass is in the 
gas; the stars contribute to 25% of the mass. The major- 
ity of the stars were actually formed in situ in the tail. 
The old stellar population which was tidally ejected along 



We computed the gas mass in a tail as the mass at radii 
larger than 30 kpc from the galaxy center. We chose to con- 
sider for the mass of the TDG progenitor all the matter in- 
cluded within a fixed aperture, rather than using a friend-to- 
friend algorithm that would better delineate the bound tidal 
object. Indeed this is how the luminosity and mass are usually 
measured in real observations. Thus, the real mass may have 
slightly been overestimated, as in the reality. The initial limit 
of 20 kpc was determined varying the aperture radius, which 
is also the method used for aperture photometry. 



i — ' — r 




Fig. 2. Bottom: contours of the gaseous component su- 
perimposed on a grayscale representation of the stellar 
component at t = 470 Myr (simulation A). Top: zoom-out 
with the contours of the dark matter haloes superimposed. 
The system is seen edge-on. 



with the gas contributes to only resp. 5% and 12% of the 
total mass. 
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Fig. 3. Contours of the gaseous component superimposed 
on a grayscale representation of the stellar component at 
t = 850 Myr (run A). The system is seen face-on. 



At t —500 Myr, the stellar disks have already merged. 
The two major condensations at the tip of the tails have 
become more compact due to the action of self-gravity. 
Smaller sub-structures developed along the two major 
tails as well as in the two smaller counter-tails. They do 
not survive later on either because they are tidally dis- 
rupted or because they simply fall back. At t =840 Myr 
(430 Myr after the periapse), the two TDGs, which now 
appear isolated, orbit around the merger remnant at a ra- 
dius of about 100 kpc. Their respective masses are now 
2.2xl0 9 M and 1.8xl0 9 M 2 . At the end of the simu- 
lations, the TDGs have converted a large fraction of their 
gas into young stars (see Fig. The mass of old stars 
generally represents 5 to 15 % of the mass of the TDG 
(see Table EJl , but this quantity may vary with the extent 
of the gaseous disk compared to the stellar one, which has 
not been changed here. 

As many real interacting systems for which detailed 
studies of TDGs exist are observed close to edge-on, we 
present in Fig. [2] an edge-on view at a time when the 
TDGs have already formed. With this line of sight, the 
tidal tails appear as linear structures and the TDGs 
as prominent overdensities at their apparent tip. This 
is close to what is observed in systems like NGC 4 676 
llHibbard fe van Gorkoml Il996ft or Arp 105 l|Duc et alJ 
Il997l) . The simulated systems seem however to have more 
narrow tidal features, especially in the stellar component. 
The fact that the mass is dominated by the gaseous com- 
ponent is also quite realistic. The HI mass of typical 
TDG candidates is indeed a factor of 2-5 higher than 



2 One fourth on the mass of the condensations of gas mea- 
sured at t —470 Myr thus turned out not to belong to the 
gravitationally linked object. 



the stellar mass es timated from their optical luminosity 
i Braine et al"]l200l|) . 

Fig. El presents snapshots 50 to 100 Myr after the peri- 
apse for runs B, C and D which were all carried out with 
the extended dark matter haloes. Despite their slightly dif- 
ferent initial conditions, orbital parameters and a spatial 
resolution which is twice lower, they all show at the tip 
of the tails prominent gaseous condensations with masses 
well above 1.0 xlO 9 M . Table [2 lists their properties. In 
run B, the longest tail gives birth to a gravitationally 
bound TDG, while the shortest tail exhibits a smooth ac- 
cumulation of matter at its tip which is not gravitationally 
bound and only transient. In run C, one of the major tidal 
accumulations fragmented into two pieces; the other one 
remains in the form of a single object. The two counter- 
tails (near y — on Fig. 0} also show accumulations of 
matter in their outer parts, but these ones do not sur- 
vive as they fall back quickly onto the parent galaxies. 
Run D, as well as run C, shows on top of the prominent 
TDG progenitor a few small stellar condensations, dis- 
tributed along the tail, with masses of the order of 10 8 M 
for the most massive ones; these gravitational clumps 
have si milar masses as those pro duce d in the simula- 
tions by IBarnes fc Hernaulstl lll992l) an d lElmegreen et alJ 
l|l993|h Like in lElmegreen et al.1 l|l993|h they were origi- 
nally formed in gaseous condensations and generally less 
than 20- 25% of their mass correspond to old stars ex- 
pelled from the parent spiral disk. A few dynamical times 
later, they mostly consist of young stars formed in situ 
from the gas clumps. Hence they may be the progenitors 
of the SSCs observed in real interacting systems. However, 
a detailed comparison of the simulated and real SSC is be- 
yond the scope of this paper. 

We carried out a few other simulations, truncating 
this time the dark matter haloes. Sub-structures form 
along the tidal tails, but no massive TDG progeni- 
tor is produced. Because these simulations are totally 
similar to those carried out by the other groups (e.g., 
IBarnes fc Hernauistlll992]) . they are not presented here. 

Therefore, including the dynamical friction in full N- 
body simulations and increasing their resolution did not 
change the basic results regarding the formation of TDGs 
initially presented in BDM03: massive gaseous clouds are 
still formed near the end of the tidal tail in this, otherwise, 
much more realistic model. 

3. Mechanism for the formation of TDGs and role 
of the dark haloes 

We have found that dark haloes should be very extended 
to enable the formation of massive TDGs. Yet, since a 
more extended halo contains more mass, one may wonder 
whether the total DM mass is more important for the for- 
mation of the TDG precursors than its actual extent. We 
checked that in BDMD03 with simulations in which the 
truncated DM halo had a mass similar to the extended 
one. We were then unable to produce any massive accu- 
mulation in the tidal tails. Thus, the DM halo extent - 
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Fig. 4. Results of the full N-body simulations of runs B (left), C (middle) and D (right) 50 to 100 Myrs after the 
periapse (see Table |2J ■ Gas contours are superimposed on the stellar component. All systems are seen face-on. 



Run 


Location 


ryi a 


Mtot 




^■young d 






(Myr) 


(M ) 


(%) 


*(%) 


A 


x>0 


60 


2.9xlO u 


5 % 


19 % 




x<0 




2.3xl0 9 


12 % 


13 % 




x>0 


430 


2.2xl0 9 


5% 


44% 




x<0 




1.8xl0 9 


14 % 


31% 


B 


x<0 


100 


2.3xlO a 


7% 


19% 


C 


x>0 


90 


2.5xlO H e 


11% 


16% 




x<0 




1.8xl0 9 


11% 


16% 


D 


y>0 


55 


l.lxlO 9 


3% 


14% 




y<0 




2.4xl0 9 


8% 


19% 



a time after the periapse. 
b total mass. 



c mass of the old stellar population from the parent's disk. 

d mass of the young stellar population born in situ. Note 
that this value is approximate as the measurement of the SFR 
is affected by the limited numerical resolution. 

e in two fragments. 

Table 2. TDGs progenitors. 



and hence the shape of the corresponding potential well 
- appears as the critical parameter, rather than the mass 
of the DM halo. Du binski, Mihos & Hernquist (1996) and 
iMihos et al.1 l|l998j) actually claimed that massive, dense, 
haloes are less able to form long-l ived tidal tails and hence 
proto-TDGs. lMihos et alJ l)l99Sj) noted however that it is 
actually "the shape and gradient of the galactic potential" , 
rather than the mass of the DM halo, that determine the 
formation and evolution of tidal tails. 

In order to decipher the various mechanisms relevant 
for the formation of TDGs and weight their relative impor- 
tance, we carried out a new set of restricted N-body sim- 
ulations. The full N-body simulations described in Sect. [3 
ensured that, using simpler models that can be more eas- 
ily analyzed, one would not disregard the main physical 
process. In this whole section, we denote as truncated a 
DM halo truncated at three times the stellar disk radius, 
and extended a DM halo extended up to ten times the 
stellar disk radius. 



3.1. A kinematical origin 

We first analyze the role of dissipation and self-gravity 
using a hybrid code: N-body for the gas, and 3-body 
for the other components. In 3-body codes, masslcss 
particles evolve in the rigid potential of two galaxies 
l|Toomre fc Toomrelll972|) . Each galaxy potential was rep- 
resented by a rigid rotation inside the central kpc, by a flat 
rotation curve for radii from 1 to 150 kpc, and a Keplerian 
potential at radii larger than 150 kpc. Such a profile corre- 
sponds to the case where the dark matter halo is extended. 

The orbit of these two rigid systems were computed 
assuming that: 

— the force exerted by a rigid system (Gl) on the other 
galaxy (G2) can be determined at the mass center of 
G2. 

— the dynamical friction can be modeled by a viscous- 
like fluid drag force prop ortional to the local density 
(e.g- JCombes et al1ll995|) . The force undergone by Gl 
is then assumed to be proportional to the density of 
G2 at the center of Gl /02(ti), and to V2 — v±, where 
Vi is the velocity of Gi. The force is then: F — kp2(r±) 
(v2 - Vi). The proportionality factor k was chosen for 
the dynamical friction time-scale at distance 50 kpc 
to be 150 Myr, which is consistent with the results of 
previous N-body simulations. This description of the 
dynamical friction may seem very approximated, but 
N-body codes with a limited number of particles, and 
an arbitrary distribution of matter, also provide only 
an approximated friction. 

We restricted our analysis to coplanar encounters, so that 
the gas particles evolve in a plane. In addition to the forces 
exerted on the two rigid galactic potential on gas particles 
in the frame of the three-body method, we used: 

— a bidimensional FFT code to account for the gas self- 
gravity 

— a sticky-particles algorithm to model the dissipative 
dynamics of the ISM 

Both codes are described in lBournaud fc Combed l)2002j) . 
We simulated galactic encounters with the same mass pa- 
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Fig. 5. Formation of the TDG progenitors in the simulations carried out with the restricted N-body + rigid haloes 
code. The evolution of the gas distribution in one of the merging galaxies is presented. The position of the the second 
galaxy is shown by a star at t =450 Myr. At earlier times, it is outside the field of view. Self-gravity and energy 
dissipation in the ISM are included in these simulations. The parameters of the simulation were: mass ratio 1:1, 8 % 
of gas, impact parameter 150 kpc, initial relative velocity 125 km.s -1 . 



rameters for dark matter and the same initial gas distribu- 
tion as in the full N-body simulations described in Sect. 2. 

Fig. presents an example of a sequence of merging 
equal mass galaxies. Even in these simpler simulations, 
massive gaseous TDG progenitors of more than 10 9 Mq 
are formed. Note that, in this run, the condensations are 
not necessarily located at the very tip of the tidal tails. At 
time t =450 Myr, the condensation at y > is at about 
2/3 of the tail extremity. This is also observed in real 
inter acting systems, such as the prototype merger NGC 
7252 l)Hibbard et al lll994j) . Its HI tidal tail extends to the 
West further than its TDG candidate, located at the tip 
of the optical tail. 

In Fig. El we show the effect of switching off the self- 
gravity and the gas dissipation, and hence the ability of 
gaseous tidal condensations to collapse. Large accumula- 
tions of gas are still present although they expectedly ap- 
pear less compact and will probably not be able to survive 
for a long time in the numerical model. This clearly shows 
that the formation of massive TDGs is fundamentally a 
kinematical process and not the result of local gravita- 
tional instabilities that would later-on grow through ac- 
cretion of surrounding material. The TDG progenitors are 



formed early on and already contain all the mass of the 
future TDG. Tidal forces within the potential well of ex- 
tended dark matter haloes account for the large accumu- 
lations of tidally expelled gas. Only then, dissipation and 
self-gravity take over and trigger the cloud concentration 
and collapse. 

One should note that, because of the kinematical ori- 
gin of the proto-TDGs, our results are not sensitive to 
the fine-tuning of the gaseous dissipation rate which is, in 
numerical simulations, often critical to structure the tidal 
tails. 

3.2. The role of the DM halo extent 

3.2.1. Response of a cold annulus to tidal forces 

In order to precisely understand the role of the DM haloes 
around the perturbing/perturbed galaxy, we analyzed how 
an annulus of cold material reacts to tidal forces. The sim- 
ulations were carried out with the same code and param- 
eters as in Sect. 13.11 Self-gravity has been disabled as we 
were only interested in studying the kinematical aspects 
of the TDG formation. 
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Fig. 6. Left: TDG gaseous progenitors at t =345 Myr for the collision shown in Fig. |3J Middle: same simulation 
suppressing the self-gravity in the gaseous component; right: same simulation suppressing self-gravity and energy 
dissipation. 



We first considered an annulus made of gas with zero 
velocity dispersion located at a constant radius Rq in the 
parent galaxy. It is perturbed by the tidal forces exerted 
by the companion galaxy. We examined how particles that 
were initially uniformly distributed along the annulus re- 
act to the perturbation 380 Myr after the beginning of the 
simulation, i.e. once the TDG progenitor is already visi- 
ble, but before its collapse (see Fig.0. As shown in Fig. [7] 
the annulus is transformed into a bent curve delineating 
the typical shape of a tidal tail. To make the analysis eas- 
ier, we subtracted the effect of the differential rotation in 
the parent galaxy. In that new frame, tidal tails appear as 
linear structures, i.e. they have a mean constant azimuth, 
instead of being curved. In Fig. |S1 we show the resulting 
plots in case of resp. extended and truncated DM haloes 
around the parent and perturbing galaxies. In both cases, 
the tail gets denser at its extremity, i.e. the particles look 
closer to each other in the 1-D space examined here. The 
amplitude of the overdensity, however, depends on the size 
of the DM haloes and is maximum when the haloes of both 
the parent galaxy and of the perturber are extended. The 
compression corresponds to the decrease of the radial dis- 
tance (dr) and more significantly to the tightening of the 
azimuthal separation (dO) between the particles. In our 
simulations, dO is three times smaller when the haloes are 
extended. 



3.2.2. Amplitude of the radial excursions 

We then considered the more realistic situation of a set of 
concentric annuli with initial radii Rq perturbed by tidal 
forces. We examined how the amplitude of the radial ex- 
cursion, i.e. the difference AR between the maximum ra- 
dial extent (reached at the azimuth of the tidal tail) and 
the initial radius (see Fig. |5J), varies with Rq. This will 
tell whether the 1-D overdensities observed for each sin- 
gle annulus overlap and hence form a 2-D condensation 
or whether they arc uniformly distributed along the tail, 




Fig. 7. Response of an annulus of uniformly distributed , 
kinematically cold, gas clouds to a tidal perturbation. 

canceling out any density contrast. Fig. El illustrates the 
two situations. In Fig. E3 we show the evolution of AR as 
a function of Rq for the truncated and extended haloes. 
In the first case, AR is a linear function of Ro, what- 
ever Rq. Thus, the separation between two annuli in the 
tail is amplified and the overdensities are stretched along 
the tails. On the other hand, for extended haloes, AR is 
nearly constant beyond a certain initial radius (15 kpc 
in the simulation shown in Fig. llOj l. Thus the distance 
between the annuli in the external regions of the parent 
galaxy is conserved and the overdensities are not diluted. 
Below this radius, the slope of the curve AR(Rq) is higher 
than in the truncated case, and hence, the matter is even 
more stretched along the tail. Both effects contribute to 
enhance the density contrast at the tip of the tail, and 
account for the accumulation of large quantities of matter 
there. In our simulations, the stretching along the radial 
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Fig. 8. Response of an annulus of uniformly distributed 
gas clouds to a tidal perturbation with DM haloes as ex- 
tended as ten stellar radii (green solid line) and with DM 
haloes truncated at three stellar radii (red dashed line). 
The x-axis corresponds to the azimuth after correction of 
the differential rotation, so that the tidal tail has a fixed 
azimuth equal to zero. The radii R have been scaled to 
the maximum extent of the tidal tail, Rmax. 
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Fig. 9. Schematic view of the effects of tidal perturbations 
to a series of concentric annuli of radius R for extended 
(left) and truncated (right) dark matter haloes. The annuli 
were initially regularly spaced out. 



direction, dAR(Ro) / dRo, is reduced by a factor of at least 
7 (i.e., the compression enhanced by this value) for the 5 
external kpc of the initial disk when the haloes are ex- 
tended. 
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Fig. 10. Amplitude of the radial excursions as a function 
of the initial radius with DM haloes as extended as ten 
stellar radii (green solid line) and with DM haloes trun- 
cated at three stellar radii (red dashed line). 

In summary, extended dark haloes cause, for a given 
annulus, a compression of matter at the extremity of the 
future tidal tail (smaller d6), and a smaller separation be- 
tween several concentric annuli (smaller dAR(Ro)/dRo). 
Thus, matter appears more concentrated both in the az- 
imuthal and radial directions than with truncated dark 
haloes. 

The physical explanation for the kinematical formation 
of a TDG progenitor lies in the shape of the tidal forces 
exerted on matter depending on its location with respect 
to the DM haloes. A mathematical description of the tidal 
field and its effect on test particles is given below. 

3.3. The shape of the tidal field 

The potential field outside of a truncated halo is a 
Keplerian potential, as arising from a point-like mass, 
while the potential field inside of an extended halo can be 
accounted for, with a good approximation, by an isother- 
mal sphere potential. One reason why a Keplerian po- 
tential and an isothermal sphere potential lead to so 
qualitatively different outcomes is likely linked to the 
properties of the tidal field of a spherical potential <j){r) 
iPekel. Devor. fc HetzronilEinl. 

Assuming the potential center lies at X = — ro, Y = 
and Z — in a cartesian (X, Y, Z) frame, we can evaluate 
the corresponding tidal field defined as the differential ac- 
celeration between a test particle located at an arbitrary 
position X, Y, Z and a test particle located at the frame 
origin. To lowest order in X/tq, Y/tq and Z/tq, the com- 
ponents of the tidal field read: 



<d 2 r2 4>)X 

-(l/r )(fi^)y 

-(l/r )(d r cf>)Z. 



(2) 
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From these equalities one can see directly that for a 
Keplerian potential (<ft oc r _1 ) one has 7y/Y = 7z /Z = 
— (l/2)'yx / X , while for an isothermal sphere potential (for 
which <j) oc logr), one has: "fy/Y = "fz/Z = —"/x/X. In 
particular, if one is interested in motions restricted to the 
(X, Y) plane, then the tidal field is divergence free for 
an isothermal sphere potential (i. e. inside a dark matter 
halo) , while this divergence is positive for a Keplerian po- 
tential (i. e. outside of a dark matter halo). In other words, 
the tidal field of a Keplerian potential tends not only to 
stretch fluid elements, but also to dilute them, whereas 
the tidal field of an isothermal sphere tends to stretch 
fluid elements while conserving their surface density. 

In order to further investigate this role of the tidal 
field shape, we have performed simple calculations of a 
cold bidimensional disk of test particles orbiting within 
its own extended halo represented by an isothermal sphere 
potential (i. e. with a flat rotation curve). It is perturbed 
by an external, time varying, tidal field, as would arise 
from a perturbing galaxy. 

The tidal field is expressed by: 

7 x = Xif(r)X 

1Y = A 2 /(r)F, (3) 

where the frame (X, Y) is centered on the disk and has 
its X-axis oriented toward the fictitious by-passing per- 
turber located at distance r(t) from the center of the 
disk. The two numerical constants Ai and A2 describe the 
shape of the tidal field and therefore correspond to its lo- 
cal properties, while the function /(r) describes how the 
tidal field intensity decays with the distance to the per- 
turber (/(r) oc r~ 2 for an isothermal sphere potential, 
while f{r) oc r~ 3 for a Keplerian potential), and is there- 
fore a global property of the tidal field. The Ai and A2 con- 
stants are actually imposed by the /(r) function, through 
the Eqs J2J). It is however instructive to artificially set A 
and /(r) separately in order to disentangle the different 
effects at work in exciting the tidal arms. 

Namely, we have run several calculations for which 
we either set Ai = 2, A2 = —1 (which corresponds to 
a Keplerian tidal field shape), or Ai = 3/2, A2 = —3/2 
(which corresponds to an isothermal sphere tide shape), 
everything else being fixed (in particular the f(r) func- 
tion). We find with this toy model that the tidal arms 
corresponding to this second choice are systematically 
more pronounced and longer than the tidal arms excited 
by a Keplerian tidal field. In a similar manner, we have 
performed calculations in which we either have f(r) oc r~ 2 
or f(r) oc r~ 3 , everything else being fixed. Although the 
tidal arms position angle might vary slightly between two 
such calculations, they display similar lengths and more 
significantly similar surface densities, which enables us to 
conclude that it is primarily the shape of the tidal field 
(i. e. the relative value of Ai and A2) that is important in 
shaping and structuring the tidal arms. Fig. II II illustrates 
this conclusion. It represents at the same date the surface 
density of a disk perturbed by a tidal field of the type 



given by Eq. with different prescriptions for (Ai,A2) 
and f(r). 

Wallin ( 1990) made a similar study of the structure 
of tidal tails, based on restricted 3-body simulations. He 
noted a density increase in the tidal tails caused by their 
twisting and the formation of orbital caustics on their out- 
side edge. Such caustics which correspond to the crossing 
of particle trajectories are also visible in Fig. ^2 Like a 
density wave, the twist creates a local overdensity which 
propagates outwards. Meanwh ile the amplit ude of the 
density enhancement decreases. IWallinl l|l990h explored a 
wide variety of initial conditions and orbital parameters, 
but did not explicitly examined how his results were af- 
fected by the shape of the tidal field, as influenced by the 
size of the DM haloes. He speculated that, when the mass 
distribution is extended, the orbital speed of the compan- 
ion decreases, and the propagation of the twist is delayed 
and slowed down. 

Our simple calculations however strongly suggest that 
it is the special shape of an extended halo tidal field (Ai = 
— A2) that is crucial to enable the formation of TDGs. 

These results agree with what we previously found: 
when the tidal interaction occurs outside truncated dark 
haloes, the divergent tidal field inhibits the formation of 
massive accumulations of matter. When the interaction 
occurs inside extended dark haloes, the tidal field can effi- 
ciently carry away from the disk a large fraction of the gas, 
while maintaining its surface density to a high value. This 
results in long, contrasted tidal tails, at the end of which 
further processes can lead to the formation of TDG-like 
objects. 

4. Discussion 

4.1. Top-down versus bottom-up scenario 

The formation mechanism of classical galaxies has long 
been a puzzle. Whether they form from the merging of 
smaller sub-structures, i.e. dwarf galaxies, - the bottom- 
up scenario -, or from the collapse of larger structures - 
the top-down scenario -, has been actively debated. The 
widely supported Cold Dark Matter model now promotes 
the hierarchical scenario where small objects form first, 
although a monolithic collapse is still regularly claimed, 
especially to account for the formation of the old ellipticals 
in clusters. 

Even if Tidal Dwarf Galaxies have a very different ori- 
gin, the discussion on their formation mechanism may, to 
a large extent, be brough t up in similar terms . The early 
numerical simula t ions b y iBarnes fc Herncmistl l)l992l) and 
lElmegreen et al.l l)l993l) showed the formation of bound 
structures all along the tidal tails, with typical masses of 
10 7 — 10 8 Mq. One may have wondered whether they were 
the progenitors of the more massive (and gaseous) TDGs 
observed in some interacting systems. Indeed they could 
in principle accrete surrounding tidal material, in partic- 
ular gas clouds. However, the uniform distribution of the 
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Fig. 11. Plots in the left column correspond to the de- 
cay of an extended halo potential (/(r) oc r~ 2 ) while 
those of the right column correspond to a Keplerian decay 
(f(r) oc r~ 3 ) of the tidal field. Similarly, the runs of the 
top row correspond to a local shape of the tidal field that 
is Keplerian (Ai = — 2X 2 ), while the runs of the bottom 
row correspond to a local shape of the tidal field corre- 
sponding to an isothermal sphere potential (Ai = — A2). 
By comparing these four plots, one can conclude that it is 
the relative value of Ai and A2 that primarily determines 
the arms morphology, i. e. the local shape of the tidal 
field. The grey levels correspond to the same values on 
the four plots, and scale with the square root of the sur- 
face density. Notice in particular the presence of a diffuse 
scattered population of test particles in the cases where 
the tidal potential has a Keplerian shape (top row). The 
sym bols in each panel re present the shape of the tidal field 
as in lDekel et all l|2003|) 



bound clumps along the tails argue against this scenario, 
at least for the massive TDGs found near the tip of the 
tails. 

Our simulations indicate that, whenever the dark mat- 
ter haloes around the parent galaxies are sufficiently ex- 
tended (their radii are at least 10 times that of the optical 
radius, or 150 kpc for a typical spiral galaxy), disk mate- 
rial coming from outside a given radius may pile up at or 
near the end of the tidal tails, instead of being stretched 
along them. Thus, a kinematical process, rather than a 
local growing instability, seems to be at the very origin 
of the progenitors of massive TDGs. Because the gas is 
more extended than the stars in the parent's disk, the 
proto-TDGs are initially mainly made of gas. Later on, 
the clouds collapse under the effect of gravity and start 
forming stars. At the end of our N-body simulations, the 



bulk of the stellar population is made of stars born in situ 
in the tail. In a few runs, the progenitor broke into two 
fragments with individual masses still above or around 
10 9 Mq. Therefore, the formation of TDGs is more con- 
sistent with a top-down scenario. 

4.2. Redefining Tidal Dwarf Galaxies 

Massive TDGs are not the only objects formed out of 
tidal debris. We enumerated in the introduction many 
other by-products of galaxy collisions, such as Super Star 
Clusters, Young Globular Clusters, Giant HII Complexes, 
etc ... Many of them are actually called TDG candidates 
in the literature. Overall, the definition of a TDG remains 
elusive. A basic consensual characteristic of any object 
defined as a T DG is that it sh ould be a gravitationally 
bound system llDuc et allEoOoh . although such a prop- 
erty is, observationally, difficult to assess ijHibbard et alJ 
2001). High resolution kinematical studies may help to 
find evi dences for a kinematical independence of tidal 
objects llDuc fc Mirabel! 1 19981 iMendes de Oliveira et all 
l200lHWeilbacher et al.ll2003HBournaud et alJl2004[h Now, 
not all self-gravitating structures in tidal tails have the 
required initial mass to form a genuine galaxy. Taking 
into account evaporation and fragmentation processes plus 
tidal disruption, a total mass as high as 10 9 M Q may be 
necessary. This is the typical mass of the giant HI accu- 
mulations observed near the tip of several long tidal tails. 
Less massive condensations may evolve, if they survive, 
into objects more similar to globular clusters. 

However, if all tidal objects were formed by the same 
physical mechanism, fixing a mass threshold for defining a 
TDG would appear arbitrary. Beside, a weakness of such 
a definition is that we do not a priori know how a proto- 
TDG evolves. Numerical simulations have just now the re- 
quired resolution to study the precise fate of tidal debris 
(see next section), while, observationally, no 'old', inde- 
pendent TDG, detached from t he parent galaxie s, has yet 
been unambiguously identified (|Duc et alJl2004l) . Their fi- 
nal mass and morpholo gy - a dwarf spheroidal, as sug- 
gested by iKrouoal l|l998j) ? - is then questionable. 

The results of our simulations and their modeling pro- 
vide a stronger theoretical support for distinguishing the 
massive gas accumulations near the tip of tidal tails and 
the other less massive condensations along them. Indeed, 
as stated in Sect. 14. ll they seem to have formed by two dif- 
ferent mechanisms. We therefore wish to define as proto- 
TDGs only the objects of the first category, i.e. those 
made of tidal material that was able to accumulate in 
the tail and later collapse, according to a top-down sce- 
nario. Indeed, this process is the only one able to form 
condensations with masses greater than 10 9 M , that are 
enough distant (> 50 kpc) from their parent g alaxies to 
avoid a rapid fall back ( Hi bbard fc Mihoall995|) or a tidal 
disruption. In contrast, the clumps along the tails have ba- 
sically a gravitational origin and formed from instabilities 
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in either the stellar fearnes k. Hernauisdll992j) or gaseous 
ijElmegreen et al.lll993|) components. 

Adopting this definition, one should admit that only 
specific, restrictive, conditions are required to form a 
TDG, in particular orbital parameters of the collision that 
favor the development of long tidal tails, and the presence 
of extended dark matter haloes around the parent galax- 
ies. Bo th criteria may actually b e simultaneously fulfilled. 
Indeed ISpringelfc White! |l999) noted that less extensive 
haloes make weaker tails. Another consequence is that 
only one proto-TDG per tail may be formed. However, 
it may later fragment (as seen in Fig. and generate 
a few extra TDGs. Determining the production rate of 
TDGs in merging systems is beyond the scope of this pa- 
per. It would require to systematically explore the param- 
eter space of the collisions (varying the orbital parameters, 
structural properties of the parent galaxies, their mass ra- 
tios, etc ...). We are currently engaged in such a long term 
program. 

Do these conditions exclude that TDGs and the less 
massive SSCs form simultaneously ? This is not necessar- 
ily the case. In the simulations shown in Fig. 0] (especially 
run D), a few faint clumps can be seen along one of the 
the tails. They look like the bound obje cts produced in the 
simulations bv lElmegreen et all 1(199.^ . although they are 
much less numerous in our simulations. Indeed, because a 
TDG already collects a large fraction of the tidal material 
(about 75%), there remains much less of it for the SSCs 
progeni tors. This would be co nsistent with the observa- 
tion bv lKnierman et all (j^QQ^) that massive and compact 
condensations do not cohabit within a single tidal tail. 

4.3. Evolution of the proto-TDGs 

In the preliminary low resolution simulations presented 
in BDM03, we noted that the TDGs survive for more 
than 1.5 Gyr and end up orbiting around the merger rem- 
nant, like any sate llite galaxy. This result is in agreement 
with the study of Iffibbard fc Mihosl l| 19951) who, study- 
ing the return of tidal material in the prototype merger 
NGC 7252, concluded that particles initially sent away at 
more than 100 kpc would not fall back within a Hubble 
time. Our new full N-body simulations have a higher res- 
olution, include feedback but were ended earlier. Still, we 
found that 0.5 Gyr after the formation of their progeni- 
tors, the two TDGs already appear isolated, orbiting at 
radial distances of about 100 kpc (See fig.|3J). Meanwhile, 
they have already used up several tens of percent of their 
gas reservoir and transformed them into young stars. This 
stellar population born in situ dominates the older com- 
ponent tidally stripped from the parent galaxies. Which 
confidence, however, can we give to these results on the 
evolution of TDGs? 

The resolution of our full N-body simulations is as a 
matter of fact still insufficient to firmly predict the sur- 
vival time of the TDGs. Each TDG progenitor is modeled 
by typically 10 5 gas particles and 10 4 star particles. This 



is just enough to model the destructive feedback effect 
(due to its internal starburst). A refined grid around the 
TDG would be required to properly take into account the 
tidal shear in the young galaxy (induced by the parent 
galaxy) and, possibly, the fragmentation of its progenitor. 
Increasing the resolution of our simulations towards form- 
ing TDGs, we should hopefully soon infer the real impact 
of these additional processes on their evolution and life 
expectancy. 

5. Conclusions 

In this paper, we presented a set of numerical simulations 
of various types aimed at studying the structure of tidal 
tails and their ability to concentrate prominent conden- 
sations likely to become Tidal Dwarf Galaxies. Our im- 
mediate objective was to check the results obtained by 
iBournaud et alJ (j2003l) with a simple numerical model: 
massive TDG progenitors may form provided that the 
dark matter haloes around their parent galaxies are very 
extended, at least ten times more than the stellar disk. 
Carrying out further tests, we reached the following con- 
clusions: 

— Gas accumulations as massive as 10 9 M are still 
formed near the end of tidal tails in N-body codes that 
include a large number of test particles, self-gravity, 
dissipation and feedback. Our initial results were there- 
fore robust. 

— Switching off the self-gravity and dissipation in sim- 
pler N-body simulations with rigid haloes does not in- 
hibit the gathering of tidal debris. This implies that 
the formation of massive tidal accumulations is fun- 
damentally due to a kinematical process: tidal forces 
make a significant fraction of the gas coming from the 
outer parent's disk to pile up near the extremity of 
the tail. Once enough material has accumulated, self- 
gravity may lead to the local collapse of the gas clouds. 

— Using toy models, we clarified the role of the dark mat- 
ter haloes around the parent galaxies. We first ana- 
lyzed how a 1-D annulus of cold gas that was origi- 
nally uniformly distributed responds to a tidal pertur- 
bation. We observed an overdensity at the tip of the 
tidal tail that is much more pronounced when the dark 
matter haloes around both the perturbed galaxy and 
the perturber are extended. We further determined the 
amplitude of the radial excursions exerted on a set of 
concentric annuli. For truncated DM haloes, the tidal 
material, and hence the individual 1-D overdensities, 
are stretched all along the tail. On the other hand, 
when the DM haloes are extended, this dilution only 
occurs below an initial radius. The density contrast at 
the tip of the tail is then enhanced. 

— Test particles expelled outside a truncated halo are 
subject to a Keplerian potential while the particles 
that will remain inside the extended halo feel an 
isothermal sphere potential. The differences in the re- 
sulting tidal fields account for the different efficiencies 
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with which large quantities of gas may be carried away. 
We found that the local shape of the tidal field plays 
the most important role. 

— The objects born in the prominent gas accumulations 
were formed according to a top-down scenario, con- 
trary to the less massive Giant HII complexes and 
Super Star (Globular) Clusters that also formed along 
tidal tails from growing local instabilities. Having dif- 
ferent physical origins, these different classes of colli- 
sional by-products have also different fates. 

— We wish to define as Tidal Dwarf Galaxies (TDGs) 
only those objects having a progenitor exceeding 
10 9 M Q , formed out of collapsing gaseous material that 
accumulated far enough from the parent galaxies, usu- 
ally (but not necessarily) near the tip of the tidal tails. 
Indeed, they are the most likely to survive and be- 
come genuine satellite galaxies. Such TDGs can only 
form under specific conditions, in particular collisions 
that favor the development of long tidal tails within 
extended dark matter haloes. They were observed so 
far in a limited number of interacting systems, but a 
systematic census of such objects is still to be done. 

More extensive simulations of galaxy collisions are re- 
quired to assess the cosmological importance of TDGs. 
The parameter space of the orbits and initial physical con- 
ditions in the parent galaxies should be more thoroughly 
explored to estimate the production rate of the proto- 
TDGs. Using a more refine grid would help to study the 
evolution of the TDGs and determine their life expectancy. 
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